Density of states, screening and the width of the quantum Hall plateaus 
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We present a consistent treatment of the quantum Hall effect within the electrostatic approximation. 
We derive the form of the density of states (DOS) which differs from the usual gaussian shape valid 
for non-interacting electrons. Below a crossover temperature, this DOS gradually transforms into the 
non-interacting DOS near the Fermi energy. We obtain an estimate of this crossover temperature, 
which determines the range of validity of this framework, and allows a reconciliation of the width and 
temperature dependence of the quantum Hall plateaus with experiments. The DOS dramatically 
enhances the electron-lattice relaxation rate, thus reducing the possibility of electron overheating. 



PACS numbers: 73.40.Hm,73.50.-h, 73.50.Jt 



I. INTRODUCTION 

The quantum Hall effect (QHE) is usually thought of 
as a manifestation of the absence of scattering of elec- 
trons carrying the Hall current. This is caused by the 
Fermi level falling between two mobility edges. Each 
mobility edge in an infinite disordered system is the en- 
ergy of a unique extended state on a given Landau level 
(LL). It has beenpBlipwn, numerically and in a number 
of approximationaJTa, that the density of states (DOS) 
of non-interacting electrons is symmetric around the ex- 
tended state energy of each LL, and has a width that 
depends on the strength of disorder but which is al- 
ways much smaller than the cyclotron energy Hlo c , with 
lo c = eB/mc. This density of states proved to be consis- 
tent with, the qualitative explanation of the widths of QH 
plateausu and with the numerical values of critical expo- 
nents describing the scaling at phase transitions between 
QH plateauaJ. Interaction between electrons has, how- 
ever, been left out of those treatments, without a strong 
a priori justification. _ 

On the other hand, EfrosU has suggested that, in pres- 
ence of long-range fluctuations of the disorder poten- 
tial, interactions between electrons are crucial, leading 
to collective screening of the disorder potential. Efros 
showed that, if the electrostatic approximation is appro- 
priate (which is true when the screening radius is much 
smaller than any other relevant length), there could be 
two types of screening which he denoted "linear" and 
"non-linear" . Linear screening is essentially metallic- 
like: in a region with linear screening the energies of 
two-dimensional electrons are confined to an interval of 
the order of temperature around the chemical potential 
and states on the top LL are partially filled. On the con- 
trary, non-linear screening is much less effective. In such 
regions the amplitude of fluctuations of the disorder po- 
tential is reduced due to electric charge of electrons, but 
this reduction is limited by the maximum density of elec- 



trons on the LL and by the inter LL energy gap. 

The existence of linear screening is supported by the 
theory of compressible edge strips in ja gate-confined two- 
dimensional electron system (2DESp. Although the na- 
ture of the potential being screened in these long com- 
pressible strips is quite different from that of the disorder 
potential, the difference in the description of screening in 
these two cases is relatively small. Conditions for the 
validity of the electrostatic approximation itself, on the 
other hand, should be universal. If such conditions are 
satisfied, one can expect a considerable fraction of the 
area of the 2DES to be covered by compressible regions 
and, therefore, a large density of states near the Fermi 
level (and at energies differing by Tiu> c ). Such DOS (we 
will refer to it as the screened DOS) is completely dif- 
ferent from the corresponding DOS for non-interacting 
electrons. The difference leads to appreciable differences 
in the predicted widths of QH plateauscl, apparently, for 
the rates of various scattering processes. Similar consid- 
erations were also successfully applied in the past to a 
quantitatively correct theory of the current breakdown 
of the QHEQ 

Dissipative transport properties at sufficiently high 
temperatures are due to thermal activation of electrons 
to unfilled LL (or due lp holes in filled LL) . It was shown 
both epqpeximentallyu'Ll and in a number of theoretical 
modelsijilj that the average dissipative conductivity of 
the 2DES in the quantum Hall regime satisfies the equa- 
tion 



Cxx = o- exp(-e a /fcT e ) 



(1) 



where ctq is generally a non-universal constant ranging 
between 0.5e 2 /h and 5e 2 //i, e a is the activation energy 
defined as the distance between the Fermi level and the 
energy of the extended state, and T e is the temperature 
of electron gas. We keep the subscript in T e to distin- 
guish it from the lattice temperature Tf, in cases when 
these are different. The range of validity of Eq. (|l|) will 
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be discussed later in this paper. We want to stress at 
this point, however, that Eq. (|l]) implies an exponential 
smallness of o~ xx , and therefore, exponentially small Joule 
heating due to scattering inside the system. 

One of the consequences of the non-interacting DOS is 
the suppression of scattering events by delocalized states 
since they are all deep under the Fermi level. Therefore, 
the energy relaxation rate via phonon emission would 
also be small in this case. If this non-interacting approx- 
imation was valid, under some conditions the phonon- 
mediated energy relaxation could become slower that the 
Joule Jj^atiELg;, resulting in observable overheating of the 
2DESl§y'EJ. Alternatively, the interacting DOS pro- 
vides for a macroscopic fraction of states at the Fermi 
level and, therefore, allows a very effective coupling of 
the 2DES to the lattice. Overheating, in this case, be- 
comes negligible. 

In this paper we address the question of the range of 
validity of both non-interacting DOS and screened DOS. 
The former is valid at low temperatures, while the lat- 
ter dominates in the higher temperature regimes. We 
present an estimate of the crossover temperature falling 
in the range between 10 mK and IK depending on vari- 
ous parameters of the sample and on experimental quan- 
tities. At temperatures above this crossover value, the 
screened DOS leads to a gigantic enhancement of the typ- 
ical phonon emission rates and, therefore, to the absence 
of overheating in the QH regime. 

The paper is organized as follows. In Sec. || we discuss 
the screening of a long-range disorder potential, closely 
following the formulation by Efroscl. The concept of 
screening is then used to introduce the new form of the 
DOS of interacting electrons. In Sec. [II we examine 
the range of validity of the electrostatic approximation 
and present an intuitive picture of the physics beyond 
this range. Section IV deals with the evolution of the 



quantum Hall effect with temperature. We derive an es- 
timate for the crossover temperature below which the 
DOS behaves more closely to the conventional DOS for 
non-interacting electrons in a random potential, at least 
close to the Fermi energy. We also propose a theory of the 
melting of the QH liquid and present numerical results 
supporting our understanding of the development of the 
DOS with the temperature. Finally, in Sec. [V] we discuss 
dissipative processes in the context of the screened DOS, 
which dramatically enhances the electron-lattice relax- 
ation rate thus eliminating electron overheating in the 
QH effect. 



II. SCREENING 

We will describe the QH system in a hierarchy of ap- 
proximations. In each hierarchical level the common goal 
will be the determination of the spatial distribution of fill- 
ing factors for each LL [n^^r)]. In a strong magnetic, 
for sufficiently slowly varying external potentials (of any 



nature) it is a good approximation to consider all the 
states on consecutive LLs to be locally separated by Tiuj c . 

The hierarchy starts with a simple electrostatic ap- 
proximation. A local filling factor is used to obtain the 
classical charge density distribution p(r) using 



2ttZ 



(2) 



N 



where I = (hc/eB) 1 / 2 ~ 100 A is the magnetic length. 
This approximation is equivalent to finding the electron 
density distribution p(r) which minimizes the free energy 
of the system. At zero temperature this corresponds to 
the minimization of a Hamiltonian consisting of the po- 
tential energy of interactions between electrons, interac- 
tion with the background disorder and quenched kinetic 
energy: 

Ul = YejS AdV JT^T + 6 / d2rp{r)V{v) 

+ J2( N +\)^c /dW^Cr) . (3) 

N - J 

The third term in the right-hand side of Eq. (|^) re- 
flects the existence of a gap between the Landau levels. 
The problem corresponding to Hamiltonian Tii was first 
analyzedi-ajidj then solved numerically by Efros and co- 
workera3ll3'E3 for a disorder potential created by remote 
donors in a plane parallel to the 2DES. As these au- 
thors pointed out, for essentially any disorder potential 
the response of the electron system results in two types 
of screening which they labeled linear and non-linear. In 
the linear screening regions the total electrostatic poten- 
tial 



$(r) 



p(v>) 



r — r 



eV(r) 



is flat and satisfies the equation 



$(r)= ti-(N+^)huj c , 



(4) 



(5) 



where /i is the chemical potential and N is the top 
partially-filled LL in the vicinity of point r. The filling 
factor n( N '(r) in such areas, which are called compress- 
ible, varies between zero and one. 

In contrast to this situation, there are also incompress- 
ible regions of non-linear screening which are character- 
ized by a complete filling of the top-most Landau level 
n^ N \r) — 1. In these cases the potential satisfies the 
inequalities 

»-(N+~)tiL> c < $(r) < n-(N+^)hu c . (6) 

Before going up in hierarchy, we give our interpretation 
of these equations. With respect to the effectiveness of 
screening, one can imagine two situations: a compressible 
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region percolates or an incompressible region percolates. 
In the former situation, there are extended states at the 
Fermi energy and the system behaves essentially as a 
metal. The latter scenario describes the quantum Hall 
effect since in this case, the extended states are away 
from the Fermi level, scattering processes are suppressed 
and the Hall conductance is quantized with a very high 
accuracy. 

A typical picture of an inhomogeneously filled 2DES is 
shown on Fig. [I]. It is the result of our numerical sim- 
ulations of the 2DES at finite temperature, the essence 
of which will be described in the following Sections. In 
the QH regime (Fig. |I](a)) the majority of the plane is 
taken by incompressible (completely filled LLs) areas. 
Yet, compressible lakes still take a finite fraction of the 
sample. By the virtue of Eq. (||) we conclude that a fi- 
nite fraction of the total number of states have energy 
equal to the chemical potential. Since all the partially 
filled regions have sizes smaller than the sample size, the 
corresponding quantum states ought to be localized and 
do not contribute to dissipation. To a limited extent, the 
situation is similar to a common 3-dimensional insulator 
system, where there are states present at the Fermi level 
but all of them are localized in space. However, there 
is a major departure from this picture: while in a com- 
mon insulator the density of states is smooth around the 
Fermi level (or even vanishes in semiconductors due to 
the Coulomb gap), in the case under study it is strongly 
peaked around the Fermi level. Finite temperature and 
neglected correlations, the effects which we consider be- 
low, tend to smear the peak, but the weight of partially 
filled states remains finite. 

Quite remarkably, the form of the DOS does not 
change dramatically with the onset of the dissipa- 
tive regime with a percolating compressible region [see 
|(a), (b) and Figs. |(a),(b)]. While the landscapes of fill- 
ing factor and potential do undergo major changes with 
the density of electrons passing through the critical value 
corresponding to the boundary between the two percolat- 
ing regimes, the DOS changes only quantitatively and in 
a smooth fashion. The existence, in both regimes, of a fi- 
nite fraction of states on the Fermi level, allows a percola- 
tion transition to occur without a visible transformation 
of the DOS. During the transition, most of the states in 
the percolating compressible region remain localized. The 
presence of these localized states at the Fermi level, how- 
ever, does have profound effects since these states can 
actively participate in low energy scattering processes, 
giving rise to finite dissipative conductivity. 

Comparing the situation above to the metal-insulator 
transition in 3-dimensional systems, we note that the lat- 
ter one occurs when the states on the Fermi level trans- 
form from localized to extended: the number of such 
states in the localized regime is so small that they are al- 
ways spatially separated from each other by wide poten- 
tial barriers and can only conduct via exponentially small 
hopping. One can say that in conventional 3D metal- 
insulator transition the Fermi level moves with respect 



FIG. 1. Electron density distribution in the 2DES; shaded 
areas represent compressible, partially filled states on the first 
or second LL, white regions are incompressible, a.: Incom- 
pressible region percolates (QH regime); b.: Compressible re- 
gion percolates. 

to the DOS; in the 2D system under consideration both 
the Fermi energy and the DOS arc almost rigid. Instead, 
the main transformation is the motion of the energy of 
the extended state, which approaches the Fermi energy 
as the system moves away from the QH regime. 

We conclude that the structure of the DOS near the 
Fermi level, in a regime when the electrostatic approx- 
imation is justified, plays a crucial role in the mecha- 
nisms of conductivity. The origin of this DOS can easily 
be seen from the Hamiltonian of Eq. [| Compared to a 
typical 3-dimensional electron system, its kinetic energy 
term is quenched by the strong magnetic field and re- 
placed by the quantized energy of electrons on a given 
LL. As long as this replacement is valid and the number 
of filled Landau levels is sufficiently small, there are no 
further degrees of freedom for the electron cloud other 
than the potential energy. In 3D systems the states with 
the same average potential energy (in the same spatial re- 
gion) would have different average kinetic energies thus 
providing for a continuous spectrum of the energies of 
these quantum states; no preference in the DOS would 
then be given to the Fermi energy. In the 2DES in a 
strong magnetic field, the potential energy is minimized 
and all the states in the "screened" region tend to have 
the same total energy, equal to the Fermi energy. 
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FIG. 2. Density of states for interacting electrons at two 
different values of the magnetic field, a.: In the quantum 
Hall regime; b.: In transitional regime. -Eo is the center of 
LL in the corresponding non-interacting system. Note that 
the only visible difference is the position of the energy of the 
extended states relative to the Fermi level. The DOS is always 
peaked at the Fermi level. 



III. RANGE OF VALIDITY OF THE 
ELECTROSTATIC APPROXIMATION 

Let us analyze more carefully the assumptions of the 
electrostatic approximation. We begin by estimating the 
size of the neglected kinetic energy term. Consider a 
mean-field equation for a single electron in the pres- 
ence of the random potential V(r) . Under the assump- 
tion of a sufficiently weak and slow varying electric field 
£ <C TiLOc/el one can expand locally V(r) and write the 
corresponding one-particle Hamiltonian keeping only the 
linear term: 



g= (p + «AAQ a 

2m 



(7) 



where £~uo/e\ is a typical local electric field due to a 
disorder potential of size uq and correlation length A. 
It is natural to choose the direction perpendicular to 
the local electric field as a quasi-translationally invariant 
one. Then, within the region where this approximation 



is valid, it is possible to enumerate the solutions by the 
wave-vectors corresponding to the quasi-translationally 
invariant-direction, as it is done in a common Landau 
equatiorO. 

This electric field produces a quadratic energy shift 
which can be referred to as the kinetic energy, since it is 
quadratic in the drift velocity £ /B. A simple estimate 
gives the value (e£l) 2 /Tilo c ^ 0.01 Tiuu c for long wavelength 
disorder potential. This justifies neglecting the kinetic 
energy term in the system we are interested in, since 
fluctuations of the bare potential commonly exceed tuv c . 
Moreover, in screened regions, the minimization of the 
potential energy leads to a simultaneous reduction of the 
kinetic term as well (decreasing of the potential fluctu- 
ation of a certain wavelength implies a decrease of the 
corresponding electric field by the same factor). There- 
fore, inclusion of the quadratic kinetic energy term does 
not produce any major change, which validates the elec- 
trostatic approximation used in Sec. [n| 

The electrostatic picture can still be changed by quan- 
tum effects. Quantum correlations contribute an addi- 
tionaL-terai into the energy density of the homogeneous 
2DESli3ta. One possible form of writing this term isE-3: 



27T 



3/2 



(8) 



where no is the density on each Landau level, v is the 
filling factor and g(v) is some complicated function of v 
of order of unity at v — 1 . This term is responsible for a 
finite screening radius r Sl contrary to the classical value 
r s = 0. It is clear that there are only two characteristic 
length scales where quantum effects set in: the size of 
the wavefunction and the average distance between par- 
ticles on each Landau level. The first length is entirely 
determined by the magnitude of magnetic field while the 
second one depends also on the density. When the filling 
factor is close to v = 1 both length scales are of the order 
of magnetic length; the screening length should, there- 
fore, be of the order of the magnetic length as well. This 
agrees with an estimate for q s = 2-n:/r s from Eq. (p|): 
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(9) 



Although Eq. (g) may not be a good approximation 
for such an inhomogeneous system as the one we study, 
it gives a good order-of-magnitude estimate, at least for 
external potentials with wavelengths \ >l. The meaning 
of q s here is the same as in the usual 3-dimensional theory 
of screening: it measures how effectively a given harmonic 
V q of an external potential is reduced due to electron 
interactions. In the 2DES with Coulomb interaction in all 
the Fourier components of external potential one has to 
make a replacement q — > q + q s in order to reproduce the 
corresponding harmonic of the screened potential. The 
Fourier component of the screened potential is F, q ~ (q + 
Qs)" 1 , and therefore screening is effective for harmonics 
with q<$^q s and much less so for q>q s - 
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The effect of electrostatic screening is the result of a 
density response while the correlation energy is deter- 
mined by the filling factor of quantum states. When 
the density and the filling factor are equal to each other, 
that is, when the relevant length scales are much larger 
than the extent of the wavefunction, electrostatic and 
correlation energies behave coherently. In regions with 
linear screening, the one-particle quantum states that 
correspond to the minimum total energy of the system, 
have energies equal to the chemical potential. Although 
electrostatic screening is limited by the effects of corre- 
lations, and therefore, the electrostatic potential differs 
significantly from the chemical potential, individual en- 
ergies, being the combinations of electrostatic and cor- 
relation contributions, fall around the chemical potential 
(ifc^T). Such states are filled partially and behave like 
a metal. 

We would like to point out that even in the regions of 
partial filling there is still a residual electric field due to 
the finite screening length. For disorder potentials with 
q > q s the electrostatic interaction becomes ineffective. 
The range of wavevectors q where, in turn, the corre- 
lation potential becomes unable to alter the individual 
particle energies prescribed by electrostatic potential is 
defined by the condition that there be one particle within 
one fluctuation of the potential. We come to the conclu- 
sion that both electrostatic interaction and correlation 
potential do not affect individual particle energies for ex- 
ternal harmonic potentials with the wavelengths A < I. 
Rephrasing it, if the system "knows" about the presence 
of fluctuations of potential with A< / , the states in this 
system are either completely filled or empty; no partial 
filling is possible. This sets a natural limit to electro- 
static approximation: if the external potential fluctua- 
tions with the wavelength A < / are important, the full 
quantum problem needs to be solved. 



IV. FINITE TEMPERATURE EFFECTS. 
GENESIS OF THE QHE 

We now consider how the picture presented in the pre- 
vious Section manifests itself in the QHE. Our object 
of attention will be so-called clean systems which are ex- 
perimentally realized in GaAs/AlGaAs heterojunctionsB. 
In those systems, a random potential is mainly created 
by remote Silicon donors placed in a layer parallel to 
the 2DES, a spacer thickness d (of the order of several 
hundred A) away from it. As demonstrated by EfrosO, 
a random distribution C (r) of such donors produces an 
electrostatic potential $(r) in the plane of 2DES which 
is given by: 

*(r) = 2tt— / ^ C(q) e l ^- qd , (10) 

e J q 

where 



Equation (|l0|) clearly shows that harmonics of the ran- 
dom potential with qd>l are exponentially suppressed. 
This lead the authors of Refs. p|, p^|Jl^ to the conclusion 
that short-wavelength potential fluctuations are not im- 
portant. Let us analyze the picture of the QHE under 
this assumption. 

For filling factors v close to an integer, say, from the 
upper side, i> = N J r 8 1 there are very few electrons above 
the top-most filled Landau level. As shown in Refs. [s|Jl5|, 
the transition between the wavelengths of the harmon- 
ics that are effectively screened (linear screening regime) 
and those which are poorly screened (non-linear screen- 
ing regime) can be written as 

R c = nf^, (12) 
o 

where C is an average areal donor density. For sufficiently 
small S , R c becomes large and surpasses the spacer thick- 
ness d. In this case some strong harmonics in the random 
potential are left unscreened since screening would re- 
quire many more electrons from the lower Landau levels 
to participate and, therefore, many excitations over the 
magnetic gap huj c . If the region with non-linear screening 
percolates, the QHE is observed since in the non-linearly 
screened regions only completely filled or empty Landau 
levels are present. For spin-degenerate Landau levels a 
rough estimate gives |#|<<5 C , with 5 C ~0.1 for the width 
of the QH plateau under typical conditions. This clearly 
contradicts with the observed QH plateaus at small tem- 
peratures which are very wide. 

Another weakness of the theory, originating in the ne- 
glect of the short wavelength harmonics of the disorder 
potential, is its prediction for the temperature depen- 
dence of the QH plateaus widthaJ. Contrary to the ex- 
perimentally observed monotonous increase of the width 
with lowering the temperatures with a feature at a non- 
universal temperature of about T ~ 300mK, this theory 
predicts narrowing of the plateaus with cooling down. 

Based on the discussion following Eq. (^) we will show 
that it is possible to complement the present theory to 
build a consistent description of the QH effect spanning 
a wide range of temperatures. Realistic conditions for 
observing the QHE assume some finite temperature T. 
This, of course, affects the density of states. The sharp 
peak at Fermi energy is smeared into a broader line. Con- 
sider a filling factor at which the region of linear screening 
percolates, |<5|>£ c - Then the states inside the percolat- 
ing region correspond to states inside the DOS peak near 
the chemical potential \i, and their energies fall within 
the temperature T around fi. This description is correct 
as long as linear screening is in effect. We have argued 
that all harmonics are linearly screened (in the sense 
explained above) as long as the wavelength X ^S> I. In 
other words, even inside a linearly screened region, fluc- 
tuations of the disorder potential with sufficiently short 
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FIG. 3. Temperature dependence of the QH plateau width. 
Measurements for twoi-samples with mobilities \x\ < fi2 are 
shown (after Gee et alH3). 



wavelengths are not significantly affected by electron in- 
teractions. However, if the amplitude of these fluctua- 
tions is much smaller than the temperature they are not 
seen by the system; they appear like small ripples on top 
of big waves [Fig. f|(a)] and the occupation numbers of 
the states involved are practically unchanged. 

When the size of these unscreened fluctuations be- 
comes comparable with temperature, however, the dis- 
tribution of occupation numbers changes dramatically, 
and eventually, when the amplitude of high frequency 
harmonics exceeds the temperature, the possibility of 
partial filling is eliminated. In this case the linearly- 
screened, metallic regions, break into pieces of completely 
filled or empty top-most Landau level. In each and all 
of these parts, the electrons behave as an incompressible 
liquid and depending on which filling factor percolates, 
the Hall conductance takes the corresponding quantized 
value [Ne 2 /h or (N-l)e 2 /h]. 

It is now possible to make a rough estimate for the 
crossover temperature. If we denote the crossover wave- 
length by A c and use Eq. ([To]) , we find that the amplitude 
V c of the harmonic A c is given by 



V c (q = 2ir/\ c ) = Voe- 2 ^^ 



(13) 



where Vq is a typical amplitude of the long-wavelength 
fluctuations. Because of the uncertainty in A c and the 
exponential dependence of V c on it, the evaluation of V c 
can be unreliable. We note, however, that for A c =d, V c = 
2 x 10 -3 Vo. Since Vo is usually of the order of several hoj c , 
and for typical densities v — 1 corresponds to Tilo c ~ 100K, 
we have an estimate of several hundred millikclvin for V c . 
As discussed above, V c has also a meaning of crossover 
temperature T c , and we therefore conclude that as we 
lower the temperature around the value T~T C given by 
the right-hand side of Eq. (|l3|), the system crosses over 
from dissipative to non-dissipative, quantum Hall regime. 
This cross-over temperature is clearly non-universal since 
it depends on the magnetic field, the spacer thickness and 
the strength of disorder. 




FIG. 4. Fluctuations of the self-consistent electrostatic po- 
tential in regions of linear screening. Long wavelength fluctu- 
ations are screened to a size of the order of the temperature. 
Short wavelength fluctuations are not screened, a.: T > T c , 
short wavelength fluctuations are negligible compared with 
T; b.: T < T c , temperature is small compared with the un- 
screened fluctuations. 



We will show below that above the crossover temper- 
ature the width of the QH plateau is almost indepen- 
dent of the temperature with a slow growth tendency 
towards lower temperatures. It is clear that upon cross- 
ing T c down the width of the plateau starts growing fast 
since more short-wavelength harmonics become "seen" 
by the system and larger areas of previously metallic re- 
gions break into incompressible pieces. Such a crossover 
feature has been-seenJn several experiments at temper- 
atures below 1KOO'E£I. At even lower temperatures the 
long-range part of disorder potential becomes completely 
irrelevant and one can treat the system as composed of 
independent electrons in a short-range disorder potential. 
It is believed that at zero temperature such a system-ex- 
hibits sharp transitions between the QH plateausE3 so 
that there is no percolating metallic region at any fill- 
ing factor. For sufficiently narrow transitional region be- 
tween the plateaus the behavior becomes universal, i.e. 
independent of the particular disorder configuration and 
electron density. This feature is entirely the consequence 
of the short-range nature of the random potential. There- 
fore, in order to approach the universal regime the tem- 
perature needs to be reduced at least below the crossover 
temperature. For very clean samples Jlor example, with 
a spacer thickness of around 1000 AEa) T c can become 
less than lOmK and universal behavior is usually not 
achieved. 

At this point, it is important to emphasize that, for 
the relatively clean samples under consideration, the QH 
plateaus are wide only at temperatures T ~ T c , where 
the quantum effects mentioned above are dominant. In 
other words, for a wide range of magnetic fields, the QH 
effect is only observed in the low temperature regime, 
and as we increase the temperature above T c , a metallic 
region percolates, extended states are close to the Fermi 
energy and dissipative processes become important. 

Consider now the situation when the filling factor is 
outside the plateau at T > T c . The extended state is then 
in the compressible region (within Ep ± ksT), and all 
states in the incompressible region are localized. Let us 
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now consider the situation in the low-temperature regime 
(T < T c ). Since states in the formerly incompressible 
region are not majorly affected by quantum effects, we 
must conclude that any extended state must come from 
one of the formerly partially filled states. It immediately 
follows that the distance between the energy of extended 
states and the Fermi level, or the width of the conven- 
tional gaussian DOS in this "quantum regime" cannot 
exceed kBT c . 

There is one additional dramatic consequence of the 
constraints on the QH liquid in the "quantum regime" : 
there is a maximum "wavelength" of a few magnetic 
lengths for an alternating structure of completely filled 
and completely empty states. Consider the transforma- 
tion of the compressible liquid as the temperature is low- 
ered past T c . Only states in the compressible regions 
can vary the occupation numbers, but these variations 
are constrained by the necessity to keep Coulomb energy 
close to the minimum. In order to achieve this, the occu- 
pation numbers must alternate with a period comparable 
to the magnetic length, so that the charge density is not 
changed significantly. This is in a deep contrast with a 
conventional perception of the incompressible QH liquid 
as homogeneous on almost macroscopic length scales. 

Similar arguments can be used to study the fractional 
QHE (FQHE). In general, one would have to talk in terms 
of composite fermions rather than electrons. This would 
complicate the picture: the Hamiltonian for composite 
fermions contains a term which makes the effective mag- 
netic field dependent on the local charge density. How- 
ever, in the incompressible liquid with a quantized filling 
factor (i.e. v = 1/3), the density is constant. In incom- 
pressible regions the disorder potential is unscreened and, 
therefore, the electric field in these regions is of the order 
of the bare electric field (~ 0.1V//um). Since the FQHE 
gap for the strongest fraction does not exceed lmeV, the 
width of the FQHE liquid may not exceed 2-3 Ib- This, 
most probably, makes it impossible to observe narrow 
electrostatic-type FQHE plateaus above the correspond- 
ing crossover temperature. It also creates additional re- 
quirements on a consistent theory of the FQHE. Since the 
FQHE is thought to be a collective effect, only the liq- 
uid with a certain degree of homogeneity can display the 
FQHE: the narrow rivers mentioned above do not satisfy 
this property. Only few electronic states across the river 
may simultaneously be the part of FQHE liquid. These 
speculations demonstrate the need for modifications in 
the Laughlin state to consistently describe the FQHE in 
a disordered system. 

Different approach has to be used in "dirty" samples 
realized in InAs/AlInAs heterojunctions or in Si MOS- 
FETS. Due to the different origin of the disorder poten- 
tial, the short-wavelength harmonics are not suppressed 
exponentially and universal behavior can, in general, be 
observed at any temperature where the system exhibits 
the QHE. If a QH plateau exists in such a system it is al- 
ways wide, and temperature induced onset of dissipation 
occurs at temperatures which have little dependence on 




FIG. 5. Illustration of the melting mechanism of the QH 
system. I stands for "incompressible", C for "compressible". 
Thermal fluctuations increase the energies of states near the 
compressible regions. 



the filling factor within the plateau. 



A. Melting of the QHE liquid 

This temperature dependent picture can be extended 
to the other limit of temperature r » T c . Then only 
long- wavelength fluctuations of the potential are relevant 
in both linearly and non-linearly screened regions. Since 
the QHE is observed when the latter one percolates, one 
may ask a question of what are the possible mechanisms 
that destroy the non-dissipative regime as the tempera- 
ture is increased. The amplitude of non-linearly screened 
harmonics is of the order of hto c while the root mean 
square of the electric field is the same as for bare disorder 
potential £ — 0.1fko c /el. If excitations over the magnetic 
gap were a relevant mechanism for melting, the corre- 
sponding temperature T m would be of the order of huj c 
while the QHE disappears completely at a temperature 
an order of magnitude smaller. The following mechanism 
gives the right order of magnitude estimate for T m . 

Consider two isolated metallic lakes separated by an 
incompressible region (Fig. ||). Energies of the states 
in the lakes can be considered to be equal to the chem- 
ical potential /Lt; states inside the incompressible region 
are completely occupied and have energies lower than /x. 
Consider a state in the incompressible region nearest to 
one of the lakes [Fig. ||(a)]. Its energy relative to /i 
can be estimated as e = —e£l, since the average extent 
of a state is of the order of the magnetic length I. If the 
temperature T is comparable to e, the state becomes par- 
tially occupied and thermal fluctuations force its energy 
to raise closer to the chemical potential. This can hap- 
pen only due to interactions with other electrons both in 
the lake and in the incompressible region. As a result of 
these interactions the size of the completely filled region 
shrinks and all states raise their energies as well. This 
situation repeats now with the new nearest-to-the-lake 
state [Fig. |(b)] . 
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It might happen that, during such a chain process, the 
electric field close to the boundary increases. Then the 
process is locked in this region unless the temperature 
is further increased. A natural estimate for the thermal 
breakdown temperature of the QHE is then 

T m = eSl = Q.lhuc . (14) 

This formula gives a range from a few to several degrees 
Kelvin, in clean samples, in qualitative agreement with 
most experiments™!. 



B. Plateau width dependence on the number of the 
filled Landau levels 

It is intuitively clear that the width of the QH plateau 
should be a decreasing function of the number of the 
filled Landau levels. To understand such dependence we 
notice that it is the charge density (not the filling fac- 
tors) which is responsible for screening of the external 
potential. Consider the same sample at different mag- 
netic fields: in the first case the completely filled first 
LL percolates, in the second it is the filled second LL 
(%LOc > Tilu^). The charge density distribution in the 
first case leads to the unscreened fluctuations of the po- 
tential of the order of fioj^ . However, such deep fluctua- 
tions are impossible in the second case: as soon as the po- 

(2) 

tential falls below fiuic ' , the 3rd LL starts filling. There- 
fore, the regions with partially filled 3rd LL at smaller 
magnetic fields are larger in size than the correspond- 
ing compressible regions with partially filled 2nd LL at 
stronger magnetic fields. Additional electrons on the top 
3rd LL must have come from some other regions in the 
system: overall it should be neutral on sufficiently large 
length scales. The only possible reservoir for such elec- 
trons are the percolating incompressible regions. By re- 
ducing the density, these incompressible regions decrease 
in size and part of them become metallic. We conclude 
that both the regions with partially filled 2nd and 3rd 
LL's are larger than their counterparts (1st and 2nd LL's 
respectively) at stronger magnetic fields. Therefore, the 
QH plateau becomes narrower with decreasing the mag- 
netic field, as expected. 



C. Density of states as a function of temperature 

We now return to the DOS and seek to understand 
how it changes when the temperature is varied. At high 
temperatures, T^>T Cl there is a peak around the Fermi 
energy with a width of the order of temperature (Fig. 
||). The amplitude of any unscreened harmonic in the 
linear screened regions is negligibly small, as compared to 
the temperature and cannot be seen inside the resonance 
curve. The tail in the DOS corresponds to states in the 
incompressible regions. 



For temperatures T c « T « T m this DOS is practi- 
cally temperature independent (except for a rounding of 
the peak at higher temperatures) . According to the pre- 
vious discussion, as T approaches and crosses below T m , 
more states from the tail of the DOS are pumped into 
the peak. The peak width, of the order of the tempera- 
ture, changes very little compared to the total width of 
the DOS. Therefore, in this range of temperatures, the 
DOS changes by increasing the height of the peak at the 
expense of the reduction of the weight of the tail. 

On the other side, for T <C T C) the amplitude of un- 
screened fluctuations becomes greater than the width 
the peak would have at this temperature in the absence 
of short-range disorder. The width of the peak in the 
DOS in this case is determined by the strength of short- 
wavelength disorder potential and is independent of tem- 
perature. This is a common non-interacting symmetric 
DOS with a few modifications: the existence of a resid- 
ual DOS at the tails spreading out to energies comparable 
with hto c , and the distance between the line of symme- 
try of the DOS and the Fermi level is of the order of 
T c <C Tilo c . The Fermi energy is either below or above 
the center of the peak depending on the filling factor. 

We were able to simulate numerically the state of the 
system for T >T C . The procedure consists of minimiza- 
tion of the free energy with respect to local density at 
finite temperature. The free energy of the system can be 
written as 

T = H! + k B T f d 2 ^ (™ W (r)lc.g[« (Ar) (r)] + 

J N 

+ [l-nW(r)]log[l-« (Ar) (r)]) , (15) 

where Hi is defined in Eq. (||) . The last term in Eq. jl5| ) 
is the entropy of the system described by Fermi statistics. 
Fig. H demonstrates the density distribution in clean sys- 
tem with average occupation number around v = 1.92. 
Part (a) corresponds to a temperature T = IK < T m . 
Figure |6|(b) corresponds to a slightly higher temperature, 
T = 2K: some expansion of the compressible region is ob- 
served. Finally, Fig. ||(c) represents the same system at 
T = 3K which is apparently greater than T m . We could 
bracket T m to within an interval between 2.5K and 2.8K. 
We also observed very weak dependence of T m on the 
filling factor. An interesting result produced by these 
calculations is that when the melting occurs at v > 2.0, 
a metallic percolating region is formed both on the first 
and second spin-degenerate Landau levels practically at 
the same time. All the features discussed above are well 
seen on these plots. Unfortunately, the electrostatic ap- 
proximation by its nature does not allow us to probe the 
system at temperatures around and below T c . A more 
sophisticated approach is needed to quantify and, prob- 
ably, numerically illustrate the crossover picture. At the 
moment we are unable to present a satisfactory, quanti- 
tative model for the microscopical nature of the crossover 
regime as well as microscopic picture of inhomogeneous 
QH liquid at temperatures below T c (see Refs. |l] pij| ). 
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T= IK, v = 0.96 




T = 2K, v = 0.96 




T = 3K, v = 0.96 

FIG. 6. Numerical simulations of the melting of the QH 
system. The same sample, with a filling factor v — 1.92, is 
studied at different temperatures. Shaded areas are compress- 
ible, white areas are incompressible, a. T = IK; b. T = 2K; 
c. T = 3K. As the temperature is raised compressible areas 
finally percolate, destroying the QH regime. 



V. DISSIPATIVE PROCESSES AND ENERGY 
RELAXATION IN THE QH REGIME 

The characteristics of the DOS at temperatures above 
T c has various implications on the dissipative dynamics in 
QH systems. We will concentrate on one of them, namely, 
the energy relaxation rate in the QH regime. By that we 
imply that we consider a "clean" system at T > T c and at 
such filling factors that the Hall conductance is quantized 
with high accuracy. In this Section we intentionally avoid 
calling this regime non- dissipative since it is dissipation, 
however small it might be, that we are interested in. 

Assume that current is injected into the system. Since 
dissipative conductivity is finite, there is finite Joule heat 
released into the electronic system. Depending on the 
rate of relaxation of this energy from the electrons to 
the lattice, the 2DES will either be in quasi-equilibrium 
(T e ~T ; ), or the electron temperature will raise (T e >Ti) 
to enhance this relaxation and achieve a steady state. 

The problem of overheating of the 2DES in such sys- 
tems has been addressed before in the study of differ- 
ent phenomena. Chow et alEj considered this effect in 
the transitional regime between QH plateaus, and argued 
that the effective electronic temperature dependence on 
the current contains needed information to reproduce 
one of the scaling exponents in the universal regime. 
Although the temperatures studied in Ref. |5J are ap- 
parently low (between lOOmK and 500mK) they may 
still be higher than T c since the sample they considered 
had a thick spacer (d = 1800A). The authors argued 
that the energy relaxation rate to acoustic (piezoelec- 
tric) phonon modes scales, at these low temperatures 
and strong magnetic fields, as ; it was implicitly as- 
sumed that the energy relaxation within the 2DES due to 
inter-electron interaction was sufficiently fast and (with 
some justification) that diagonal conductivity a xx was 
independent of wavevector q and electron temperature 
T e . The latter assumption obviously fails in the QH 
regime since finite conductivity results from either ther- 
mal activation (~ exp[e a /fcsT e ]) or variable range hop- 
ping (~ exp[— A/T@]), where e a is an activation energy, 
and (3 is some mo del- dependent positive power. 

Our analysis of heating at filling factors within the QH 
plateaus is reminiscent of the analysis of Ref. ^J. How- 
ever, we will use explicitly the properties of the DOS, 
since it is unclear how to express the surface acoustic 
wave absorption probability through a global conductiv- 
ity in the case of spatially inhomogeneous system (see 
also Ref. ^B|). 

The notion of overheating of the electron system as- 
sumes that the electrons are considered to be in ther- 
mal equilibrium at some temperature T e . This implies 
that all the energy pumped into the system is distributed 
throughout it at rate t~ 1 which is much greater than the 
rate of any other energy relaxation process in the system. 
Important processes which are to be considered are the 
supply of energy from some external source and release 







of this energy to a thermal bath. If r e is the shortest 
time scale in the problem, and if the temperature of the 
thermal bath is denoted as Tf,, the 2DES temperature T e 
is defined by the energy balance equation 



Pin -fout(-^e) ; 



(16) 



where Pj„ is the input power, and P ou t is the energy re- 
laxation rate from the electrons to the heat bath. The 
relevant question is how much the electron temperature 
needs to be raised so that Eq. [l6] is satisfied. If energy 
relaxation to the lattice is slow, it is clear that consider- 
able overheating (T e ^>Tf,) will be observed. Conversely, 
a fast relaxation mechanism assures that both tempera- 
tures are similar (T e ~Tf,). 

The energy fed into the 2DES QH system is Joule heat 
associated with the injected current j and dissipative re- 
sistivity p xx : Pi n — p xx j 2 . This assumes that most dis- 
sipative processes do not involve electron-phonon inter- 
actions, otherwise part of the heating is directly applied 
to the lattice. At this point it is important to emphasize 
that in the QH regime p xx is exponentially small. For 
the phonon emission power we write: 

P out (T e ) = 2tt^ Y, Mq) I M i!f (q) | 2 /, (1-//) x 
ij q 

[(n q +l)<5[e»-e/-7iw(q)] - n q <5[ej-e/+7ia;(q)]] . (17) 

where indices i and / correspond to initial and final elec- 
tronic states respectively, q is the phonon momentum; 
fij(T e ) are fermionic filling factors at temperature T e , 
n q (Tb) is a bosonic filling factor at lattice temperature Tf,, 
and Mij(q) denotes the matrix element of the electron- 
phonon interaction. In Eq. ( p"7j ) the first term inside the 
vertical brackets stands for phonon emission while the 
second corresponds to phonon absorption. To estimate 
-Pout we make the following approximation: each elec- 
tronic state can be visualized as a strip with a transverse 
size of the order of the magnetic length / stretched along 
some equipotential. Inside any individual potential fluc- 
tuation this potential can be linearized. The wave func- 
tion in this linearized region can then be written as a so- 
lution of Landau equation in a homogeneous electric field. 
Solutions of this equation arc enumerated by the local set 
of one-dimensional wave vectors k that have meaning of 
the position of the wave functions centered, which for the 
lowest Landau level read 



ipk(x,y) = 



1 



7T 1 / 4 



exp 



ikx 



(y - yo? 

2/ 2 



(18) 



where x is the direction along the trajectory and y$ ~ kl 2 
is the guiding center position. Energies of these states 
follow the external potential in the fluctuation: 



1 , 



(19) 



with S being the typical electric field inside the fluctu- 
ation. Since the characteristic wavelength of potential 



fluctuations in the regime we are studying is much larger 
than I, and any transition due to electron-phonon in- 
teraction is possible as long as initial and final states 
have substantial overlap, such transition occurs when 
the electron moves from one strip to another approxi- 
mately I away in the same potential fluctuation. Using 
the dispersion relation for longitudinal acoustic phonons 
lu = c^q, with q = [q 2 + q 2 + Q 2 ] 1 / 2 , an estimate of 
£~0.1 huj c /el, and given that cl ~ 10 3 m/s, one can make 
sure that single phonon processes are indeed allowed: 
e£Ayo = CLq>CLq x ^CLAyo/l 2 leads to £>0.005 hu> c /el, 
which is certainly the case. The matrix element Mij(q) 
calculated with the wave functions from Eq. ([l8]) depends 
on the wave vector of initial state only through a phase 
factor that is irrelevant in this problem. 

Since each strip passes through many potential fluctu- 
ations (provided that the corresponding equipotential is 
sufficiently long) the electron on the strip probes many 
random configurations of the potential and, therefore, av- 
erages over the whole ensemble. Hence, one can use the 
ensemble averaged density of states -D(e) to describe the 
probability of the electron having particular initial and 
final energies, at least for electrons in the incompressible 
regions. In this approximation, Eq. ( |l7j ) takes the form: 



P™t= 2tt 



E 



M(q)| 2 
D(e k )D{e k - 



1 



q 



a- 1 
- c L q) 



(20) 



(P + x)(l+x) 



Six 



D(e k )D(e k +c L q) 



where we denoted x — exp[e k /kT e ], a = exp[cLq/kTf,] and 
/3 = exp[cLq/kT e ]. Changing to an integral over e k and 
after simple manipulations we obtain: 



= k B T e 



|M(q)P 



a — [3 
a-1 



[°° D{k B T e \ogx)D{k B T e \ogx/P) 
Jo X (f3 + x)(l + x) ' 1 ] 

where we absorbed all numerical factors into a redefini- 
tion of the matrix element A/(q). We did not have to 
specify the exact form of electron-phonon interaction; it 
can, in fact, be of either piezoelectric or deformation po- 
tential nature. The relative contribution of these two 
types of interaction depends on the temperatureEH, how- 
ever both are of the same order of magnitude at temper- 
atures around IK. After averaging over phonon polariza- 
tion, one gets matrix elements of the order of unity. We 
chose not to go into details of the form of these matrix 
elements (see Refs. |27|-|30|) since for our purposes only the 
order of magnitude matters. We kept the factor k B T e as 
it appeared with the change of integration variable for 
the sake of dimensionality. 

The most important part defining the magnitude of 
P ou t in Eq. ( pl| ) is the density of states D(e). For exam- 
ple, if only states well below the Fermi level, with ener- 
gies e k B T e , have weight in the DOS then the range 
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of integration over x in Eq. ( pl| ) becomes exponentially 
narrow: the upper limit of integration is then of the or- 
der of exp[— eo/ksTe], where eo is the highest energy such 
that D(eo) is not exceedingly small. On the other hand, 
if the DOS spreads all the way to the Fermi level, the 
corresponding integral in Eq. (|2lj) is not small nor is the 
emitted power P ut- 

To see that energy relaxation rate due to phonon emis- 
sion is sufficient to keep the difference T e — negligi- 
bly small, we recall that Joule heating rate is suppressed 
exponentially by a factor of exp[— e a /kT e ]: dissipative 
transport is due to the electrons (or holes) in the ex- 
tended states. Different mechanisms of scattering are 
available for these electrons. Scattering off the electrons 
in metallic, compressible regions is one mechanism lead- 
ing to thermal equilibration within the 2DES. The other 
mechanism proposed intiJO involves interaction with lat- 
tice phonons. It is not clear which of these contribu- 
tions to dissipative resistivity prevails. If under certain 
conditions the latter one dominates, the 2DES has to 
be in thermal equilibrium with the lattice and T e ~ 
Electron-electron scattering is the only way of releasing 
energy into 2DES. It is clear, however, from Eq. (21), 



as well as from intuitive arguments, that electrons in 
metallic regions are interacting with the phonons very 
effectively and are never a bottleneck for thermal equili- 
bration with the lattice. In general, the metallic regions 
are always in thermal quasi-equilibrium with the lattice. 
They may not be equilibrated with the electrons in in- 
compressible regions: in this case the QH system has 
to sustain the temperature gradients. However, it as- 
sumes that the dissipation of energy by incompressible 
electrons carrying the Hall current occurs via interaction 
with phonons: in the QH system diffusion (and therefore 
dissipative conductivity) occurs through inelastic scat- 
tering due to the non-degenerate nature of the states in 
long-range disorder potential. 

Expression (|l]) is applicable in the transitional regime 
as well. However, when the metallic region percolates, 
dissipation is also greatly enhanced. Both Joule heating 
and phonon emission are determined in this case by the 
same DOS, that is by a finite fraction of the total number 
of states at the. Fermi level. Observation of overheating 
in this regimetj is therefore not surprising. 



VI. CONCLUSIONS 



In this paper we demonstrate that screening of long- 
range fluctuations of the disorder potential leads to 
a novel form of the density of states. Contrary to 
the conventional DOS for non-interacting electrons, this 
screened DOS possesses a distinctive peak at the Fermi 
level originating from states in metallic (compressible) 
regions. We argue that this form of DOS remains al- 
most intact with changes in the magnetic field, the only 
variation being in the width of the peak and in the posi- 



tion of the energy of the extended states relative to the 
Fermi level. We also argue that, as the temperature of 
the system crosses below a certain value T c , short-range 
fluctuations of the disorder potential become dominant 
and the nature of the states within the narrow peak in 
DOS changes. This T c is the lower boundary of the range 
of validity of the electrostatic approximation. 

For filling factors that only exhibit the QHE at low 
temperatures (T <~ T c ), short-range fluctuations domi- 
nate and the full quantum problem needs to be solved. 
Based on numerical and analytical results obtained by 
a number of authors, one can expect that for non- 
interacting electrons, the DOS has symmetric gaussian 
shape, centered at the energy of extended states. In light 
of our developments, the assumption of non-interacting 
electrons needs to be amended by the requirement that 
the charge density remain unchanged when the temper- 
ature is lowered through T c , that is, it is still prescribed 
by electrostatic solution even at low temperatures. This 
leads not only to a constraint on the distance between 
the energy of extended states and the Fermi level never 
to exceed T c but also to a more dramatic conclusion: 
namely, that in these, short-range fluctuations dominated 
regime, the percolating quantum Hall liquid can nowhere 
be wider than few magnetic lengths. This conclusion calls 
for a closer look at the existing quantum theory of the 
QHE which is based on the assumption of an homoge- 
neous system. 

The current common perception of the QHE is that 
it is an entirely quantum phenomenon. In this sense, 
the theory of the QHE refers to the lower temperatures 
(T < T c ) and to the magnetic fields away from the cen- 
ter of the plateaus. Yet it is important to bear in mind 
that this QH system is very inhomogeneous and that the 
"pure" considerations do not apply to the whole system 
but rather to the regions where the disorder potential 
is metallically screened. All the observed quantum ef- 
fects, including the universal behavior at low tempera- 
tures, have their origin in these regions. The complete 
theory of the IQHE has to account for this non-trivial 
microscopic structure of the 2DES. 

Based on the understanding of the QH system de- 
scribed by the electrostatic approximation, we also 
present a model of melting of QH liquid at higher tem- 
peratures. We were able to explain why the electrostatic 
solution predicts narrow quantum Hall plateaus while ex- 
perimentally observed plateaus are usually much wider, 
this difference originating in the crossover behavior de- 
scribed above. 

Finally, we applied the derived form of the DOS to the 
study of dissipative processes in QH system. We showed 
that since Joule heating is limited in the QH effect by a 
very low dissipative conductivity, and the relaxation of 
energy due to the coupling to phonons is very effective 
(because of the large number of states on the Fermi level), 
electronic overheating cannot occur (or is negligible) in 
the QH regime. 
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